Heating and Turbulence Driving by Galaxy Motions in Galaxy 

Clusters 



Woong-Tae Kim 

Department of Physics & Astronomy, FPRD, Seoul National University, Seoul 151-742, 

Republic of Korea 

wkimSastro . snu . ac . kr 
ABSTRACT 

Using three-dimensional hydrodynamic simulations, we investigate heating 
and turbulence driving in an intracluster medium (ICM) by orbital motions of 
galaxies in a galaxy cluster. We consider A^^ member galaxies on isothermal and 
isotropic orbits through an ICM typical of rich clusters. An introduction of the 
galaxies immediately produces gravitational wakes, providing perturbations that 
can potentially grow via resonant interaction with the background gas. When 
Ng^'^Mii <, 100, where Mu is each galaxy mass in units of 10^^ Mq, the pertur- 
bations are in the linear regime and the resonant excitation of gravity waves is 
efficient to generate kinetic energy in the ICM, resulting in the velocity dispersion 
(T„ ~ 2.2Ng^'^Mii km s~^. When N^^Mu ^ 100, on the other hand, nonlinear 
fluctuations of the background ICM destroy galaxy wakes and thus render res- 
onant excitation weak or absent. In this case, the kinetic energy saturates at 
the level corresponding to ~ 220 km s~^. The angle-averaged velocity power 
spectra of turbulence driven in our models have slopes in the range of —3.7 to 
—4.3. With the nonlinear saturation of resonant excitation, none of the cool- 
ing models considered are able to halt cooling catastrophe, suggesting that the 
galaxy motions alone are unlikely to solve the cooling flow problem. 

Subject headings: cooling flows — galaxies: clusters : general — turbulence — 
waves — X-rays: galaxies 



1. INTRODUCTION 



A lack of cold gas in the central parts of rich galaxy clusters, as revealed by high- 
resolution X-ray observations, has posed a "cooling flow" problem, requiring sources of heat 
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to balance radiat ive cooling of an intraclu ster medium (ICM). Among the proposed heating 
mechanisms (see iPeterson fc Fabianll2006l for review), energy injection from active galactic 
nuclei appears to be the mos t favorable, although it requires a central black hole to be 
more massive than observed (iFujita fc ReiprichI 12004') and is una ble to maintain a long 
term energy balance if jets are narrow (jVernaleo fc Reynolds 2006 ). Diffusive heating via 
thermal conduction and/or turbulent mixing may also be effective i f the relevant diffusion 



coefficient is quite large a nd fine tuned (e.g., iKim fc NarayanI l2003l : IVoigt fc FabianI 12004 
Dennis fc Chandranll2005l ). 



Less well recognized is th e ICM heating by cluster galaxies that possess a lot of ayailable 



energy in their motions (e.g., lMillerlll986l : iBregman &: David! Il989l ). iBalbus fc Sokerl (11990 



hereafter BS) showed that resonant excitations of internal waves driven by orbiting galaxies 
produce heat at cluster centers comparable to radiative los s , prov ided the galaxy mass is 
large enough. Using Monte-Carlo approaches, lEl-Zant et all (120041 ) showed that dynamical 
friction of g alaxies can be a distributed source of heat if the mass-to-light ratio of galaxies 
exceeds 10. iKim et all (120051 ) found that heating by dynamical friction reduces the growth 
rates of thermal instability. All these works suggest that the effects of galaxy motions on 
thermodynamic evolution of ICM are by no means negligible. 

Most of the studies cited above are based on the assumption that the energy lost by 
galaxies is all transferred to thermal energy of the ICM locally near the galaxies. This is 



(Deiss & Just 


1996: 


Ostriker 


1999) 



The possibility of the ICM stirring by ga laxy motions is intere sting because turbulence 



appears to be pervasive in the ICM (e.g., ISchuecker et al.l |2004| ) and p erhaps determines 



the characteristic streng ths and scales of cluster magnetic fields (e.g., IClarke et all 12001 



Subramanian et alll2006l ). Although BS allowed for the spatial propagation of inter nal waves, 
they focused on linear gaseous responses in the WKB limit. iLufkin et all (119951 . hereafter 
LBH) explored nonlinear evolution of gravity waves, but their models considered a single 
galaxy on a radial orbit. In this Letter, we extend LHB by considering a number of cluster 
galaxies on isotropic orbits. By varying the mass and number of galaxies, we quantify the 
thermal and kinetic energies induced by galaxy motions and explore the level and shape of 
such driven ICM turbulence. 



2. Model and Method 



We consider a galaxy cluster in which the ICM is initially in hydrostatic equilibrium 
under a dark matter potential $dm- For the initial temperature of the ICM, we adopt a simple 
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form T(r) = 7[1 — 0.6e (^Z^':)^] keV with the coohng radius Vc = 125 kpc. We represent a 
rigid dark halo using an NFW profile with th e characteristic mass Mq = 6 x 10^^ Mq and 



scale radius = 460 kpc (INavarro et al.l 119971 ) . Although we do not allow for the presence 



of a central dominant galaxy considered in LBH, the cooling core inside in our model still 
satisfies the condition for resonant excitation of gravity waves (see §3.11) . 

To study the responses of the ICM to member galaxies, we solve the ideal hydrodynamic 
equations: 

(5/9t + v- V)p + pV-v = 0, (1) 

p{d/dt + V ■ V)v = -VP - pV($DM + $,), (2) 

p{d/dt + V • V) (e/p) = -P V • V - A, (3) 

where $g is the time-varying gravitational potential due to the galaxies and A is the volu- 
metric cooling rate. The other symbols have their usual meanings. The effects of gaseous 
self-gravity and magnetic fields are ignored. We adopt an ideal gas law P = (7 — l)e with 
7 = 5/3. We run both adiaba tic (with A = 0) and coo li ng (w ith A 7^ 0) models by taking 



the cooling function A used in iRuszkowski fc Begelmaru (120021 ) 



We consider a total of A^^ member galaxies distributed within 1 Mpc. We assume that 
they follow isotropic and isothermal orbits with velocity dispersion ar = 800 km s~^, in which 
case the equilibrium galaxy number density is oc (1 + r/rs)^^^^"^, where r] = 2GMq/ {rga"^) 
17.5 ( Kim et a ] l2005h . Under this distribution, about 13% of the galaxies are located within 



at any given time. We ignore back reaction of the ICM to the galaxies since the energy 
lost due to dynamical friction is small (LBH). We represent each galaxy using a Plummer 
potential $p(r) = —GMg{r'^ + a^)~^/^ with mass Mg and scale radius a = 7 kpc; the total 
perturbing potential is constructed as $g(r, t) = Yl!i=i ^pd^^^^il^)!)' where is the position 
vector of the i-th galaxy at time t. To simulate diverse cluster conditions, we vary the 
number and mass of the galaxies in the ranges of A^^ ~ 10^ — 10^ and Mg ~ 10^^ — 10^^ M^^. 



]2 - 10^ and Mg ~ 10" - 10^^ 
We report in this work the results only for models where all the galaxies have equal masse; 



We follow the nonlinear evolution of the hot ICM using a modified version of the ZEUS 



code (jStone fc Norman! Il992l ). parallelized on a distributed- memory platform. Our simula- 
tion domain is a cubic box with each side of 2 Mpc; the center of the box is located at the 
cluster center. We construct a logarithmically spaced Cartesian grid with 256^ zones, with 



^By running a few mod els in which galaxy m asses vary according to the Schechter function with a mass- 
to-light ratio M/L - L°-^ |Gerhard et al.ll200lh . we have confirmed that results are almost unchanged if the 
breaking mass in the mass distribution is equal to Mg in the corresponding fixed mass case. 
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outflow conditions at all boundarie^. The grid spacing is 0.3, 4, and 33 kpc at the center, 
the cooling radius, and the edge of the box, respectively. We have confirmed the accuracy 
of the code by comparing th e test results for wakes produced by linear-trajectory perturbers 
with the analytic formula of lOstrikerl (119991 ). 



3. Nonlinear Simulations 

3.1. Resonant Excitation by a Single GalcLxy 

We first explore the responses of the adiabatic gas to a single galaxy moving either on 
a radial or a circular trajectory. BS showed that gravity waves excited by a galaxy with 
orbital frequency cUorb become trapped and amplify in the region where the local Brunt- 
Vaisala frequency u^v exceeds ci^orb- This finding was subsequently confirmed by LBH who 
ran numerical simulations for a radial-orbit galaxy. Our aim here is to find the dependency 
of energy injection rate on the galaxy mass as well as on the shape of its orbit. 

Figure [T^ shows the radial distributions of the radial orbit, circular orbit, and Brunt- 
Vaisala frequencies in our ICM model. Even without a central massive galaxy, cjbv remains 
almost flat inside 80 kpc, allowing resonant excitation of gravity waves there. For models 
with a radial orbit, the galaxy initially set to move with velocity v = 1740 km s~^ from the 
origin has an average speed of 820 km s~^ until it reaches a turning point at r = 150 kpc 
(black dot in Fig. [T^). For circular-orbit models, the orbital speed is f = 1060 km s~^ at 
r = 130 kpc (blue dot). In both cases, the orbital period of the galaxy is 0.73 Gyr, which is 
chosen to ensure cubv > cUorb inside Vc- 

As the galaxy starts to move, it generates density and velocity perturbations, forming 
a gravitational wake and imparting some of its gravitational energy to thermal and kinetic 
energies of the ICM. In general, the perturbations are a superposition of p- and g-waves. 
While high-frequency p-waves propagate out through a stratified background, outgoing g- 
waves reflected at, and remain trapped within, the resonance radius (~ Vc in our models). 
The gas near the center receives periodic kicks from the galaxy, enhancing the levels of 
density and velocity fluctuations. Overall evolution of the models with a galaxy on the 
radial orbit is similar to that presented in LBH. 

Figure [lb shows a snapshot of the perturbed density in the orbital plane at t = 6 Gyr 



^We have also run models with 128'^ zones, and checked that the results are within less than 10% of those 
from 256'^ runs. This suggests the energy dissipation caused by numerical diffusion is tolerable. See i^3.2l 
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for a model where a galaxy with mass Mn = Mg/(10^^MQ) = 5 orbits circularly at a 
near transonic speed (Mach number = 0.97) in the clockwise direction. It is apparent that 
perturbations periodically provided by the orbiting galaxy is focused to the central part. 
The associated vorticity is also well contained inside the resonance marked by a dotted 
circle, indicative of resonant excitation (LHB). The density perturbations near the galaxy 
i n a weak traihng s hape are a characteristic feature of a wake for a circular-orbit perturber 



(IKim fc Kim 1120071 ). 



As the galaxy orbits the cluster center and resonantly interacts with the background 
gas, the kinetic energy absorbed in the ICM secularly increases (approximately linearly) 
with time, while showing some temporal fluctuations. We run a number of models with 
varying Mg, and measure the rate E^^ at which kinetic energy increases inside Vc- Figure 
^ plots the resulting Ek, which are fairly well fitted by Ek = 1.2 x 10^°Mf^ erg s~^ and 
Ek = 3.5 X lO^^Mf^ erg s~^ for the radial- and circular-orbit cases, respectively. For the 
model parameters we adopt, therefore, a galaxy on a radial orbit is about three times more 
efficient in driving kinetic energy into the ICM than the circular-orbit counterpart, since it in 
the former orbit can traverse the central region directly (BS). The dependency of Ek on Mg 
indicates that perturbations in all the models with a single galaxy are in the linear regime. 



3.2. Heating by Cluster GalcLxies 

We now consider more realistic cluster models in which the ICM is continuously stirred 
by many member galaxies. We run nine adiabatic models as well as nine cooling counterparts. 
Figure [2] plots time evolution of the kinetic energy Ej. of the ICM inside for some of 
the adiabatic models. A sudden introduction of the galaxies causes the ICM to respond 
abruptly, initiating an rapid increase of Ek for t < 0.1 Gyr. For models D [Ng = 10^, 
Mil = 5) and E {Ng = 10^, Mn = 1), the perturbations are initially in the linear regime and 
soon begin to interact resonantly with the background gas. As gravity waves concentrate 
toward the center and amplify, Ek grows secularly with time. For t > 2 Gyr, Ek in models 

1 /2 

D and E (and other models with Ng' Mn < 100) increases almost linearly at a rate Ek = 
2.6 X 10^^ {Ng/ 10^) Mf^ erg s^^. Compared to the results of §3.11 this corresponds to about 
20 radial-orbit galaxies taking part in resonant excitation of the ICM inside Vc- On the 
other hand, in model A (A^^ = 10^, Mn = 10) where initial perturbations are very strong, 
Ek stays almost constant at ~ iQ^o.b during its entire evolution, indicating weak or no 
resonant excitation. Models B {Ng = 10^, Mn = 5) and C {Ng = 500, Mn = 5) that have 
weaker initial perturbations than model A enhance Ek to the saturation level ~ IQ^'^-^ erg for 
t ^ 1 Gyr, after which Ek again remains constant. At saturation, the amphtudes of density 
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fluctuations are about 40% relative to the mean value, easily disrupting galaxy wakes that 
would otherwise supply fresh perturbations for resonant excitation. Consequently, resonant 
excitation becomes weak or even absent in a highly nonlinear background. 

Unlike kinetic energy that grows with time when NgM^ is sufficiently small, we find 
that thermal energy of the ICM does not show a clear indication of resonant excitation. The 
total amount of thermal energy within driven by a mere introduction of the galaxies is 
AEt = lO^^NgMii erg, which is already large and stays more or less constant with time for 
all the models. This may be because kinetic energy is more prone to resonant excitation, 
or because the amount of heat supplied (presumably at a rate similar to Ek) by resonant 
excitation is much smaller than AEt, so that it does not readily manifest in the energy curves 
over time. The ratio of kinetic to thermal energies in a wake produced by a single galaxy 
with size a and vel ocity cr„ is roughly ~ {GMg/ a^ja)'^ / (Sc^) , where Cg is the adiabatic sound 
speed of the ICM ( Just fc Kegel 199Cll ). Since this ratio is less than 10^^ when Mu < 1 for 



typical parameters we adopt, kinetic energy is certainly a much better tracer of resonant 
excitation. 

To check if galaxy motions and the associated heating can solve the cooling flow problem, 
we repeat the simulations by including the radiative cooling explicitly. Our models lose 
thermal energy at a rate Lx ~ 10^'* erg s~^ from the cooling core. Without any heat source, 
they would experience a cooling catastrophe within 0.6 Gyr. It turned out that none of the 
models we considered were able to prevent, albeit considerably delay, the runaway cooling. 
For instance, the cooling model with Ng = 10^ and Mu = 5 (corresponding to model 
D) undergoes a catastrophic event at 2.5 Gyr, while the model A^^^ = 10^ and Mu = 10 
(corresponding to model A) develops a cooling flow in 1.7 Gyr near the center. The heating 
rate due to resonant excitation is probably lower in the latter model, for which the perturbed 
kinetic energy saturates immediately. 



3.3. Properties of ICM Turbulence 

We have seen that motions of the member galaxies generate a large amount of kinetic 
energy in the ICM. The kinetic energy is in the form of fluctuating isotropic velocity fields 
with vanishingly small mean values. We regard these spatially-uncorrelated, random gas 
motions as ICM turbulence. To quantify the turbulence level, for each run we measure 
the velocity dispersions of the gas inside Tc- Figure [3^ plots as open circles the resulting 
density- weighted 3D velocity dispersion = {2Ek/ Mc)^/"^ averaged over 3 — 6 Gyr as a 
function of NgMl^ where is the ICM mass inside Tc- Filled circles give non- weighed 
velocity dispersions, which is larger than the density-weighted values by about a factor of 
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1/2 

1.5. Consistent with the results of the previous subsection, increases hnearly with Ng Mn 
until Ng^'^Mii ^ 100, beyond which is approximately constant at ~ 210 — 230 km s~^. 

To characterize the turbulence driven by galaxy motions, we calculate Fourier power 
spectra of the compressive and shear components of ICM velocities defined by = |k ■ v^p 
and = |k X v^p, respectively. Here, is the Fourier transformed velocity and k is the 
wavenumber. We then bin them spherically in the k-space and calculate the angle-averaged 
power spectra PdK) and PsiK) as functions of the radial wavenumber kr- Figure [3t shows 
Pc and Ps for models C and E at t = 6 Gyr. The ratio of total power in the compressive 
to shearing parts is abou t 2.5 for both models, indicative of fairly subsonic turbulence (e.g.. 



Vestuto fc Ostrikeiil2003l ). The power index of the shearing part is —3.7 in the inertial range 
for model E, which becomes steeper with increasing cr^, yielding —4.3 for model d. The 
compressive parts have similar slopes, although they have excess power at ~ 50 — 100 kpc 
scales, which appear to be associated with the mean galaxy separation^ 



4. Discussion 

Galaxies in a cluster contain a plentiful amount of energy in their orbital motions 
that can be trans ferred to turbulent kinetic and thermal energies of the ICM (e.g., BS; 



Deiss fc Justlll996l ). In this Letter, we have shown that resonant excitation of gravity waves 
driven by galaxy motions is efficient only when density and velocity fluctuations in the back- 
ground are in the linear regime, while becoming inefficient when the background medium 
exhibits large amplitude fluctuations. Although it is uncertain at what rate resonant ex- 
citation heats the ICM, our numerical results suggest that heating by galaxy motions is 
insufficient to quench the cooling catastrophe. If the heating rate is similar to the kinetic 
energy injection rate Ek we found, the energy balance between the heating and X-ray cooling 

1/2 

requires Ng Mu ~ 200. Although rich clusters may contain enough number and mass of 
galaxies to satisfy this condition, resonant excitation will switch off when this condition is 
met. We thus conclude that heating by galaxy motions alone cannot be the main solution 
to t he cooling flow problem, although it can delay the catastrophic event signiflcantly (e.g.. 



BS: lKim et alll2005h 



The two key parameters that control the rate of energy injection and the level of turbu- 



Caution should be made in interpreting the power indices in comparison with Kolmogorov power spectra, 
since in our models the background density is not uniform and the velocity field is not periodic. 

*The galaxies are concentrated more strongly toward the center and have a mean distance of 82 kpc inside 
the cooling radius. 
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lence are the mass and number of galaxies. Many uncertainties surround the observational 
determinations of the average galaxy mass, but a recent analysis using strong-lensing models 
shows that a cluster galaxy can have mass as large as 5 x 10^^ Mq including a dark halo 
( iHalkola et al.l 120071 ). For rich clusters with A^^ > 10^ inside 1 Mpc, therefore, our numer- 
ical results suggest that ICM turbulence driven solely by galaxy motions is probably in a 
saturated state with a„ ~ 220 km s~^. Given the many arbitrary choices for the cluster 
pa rameters, this value is in rough agreement with an analytic estimate of ~ 300 km s~^ 
by ISubramanian et all ( l2006l ) who used scaling laws of hydrodynamic wakes without con- 
sidering nonlinear saturation. Note that the saturated is similar to those required to 
explain the observed magn etic field strength (~ lyuG) in terms of energy equipartition (e.g., 
Goldman fc RephaelilEoQll ). 



The angle-averaged velocity power spectra of turbulence driven in our models are char- 
acterized by inertial-range slopes ranging from —3.7 to —4.3. It is interesting to note that 
these are comparable to the values between —11/ 3 and —13/3 inferred from the observed 



pressure maps of the ICM in the Coma cluster (ISchuecker et al. 



inhomogeneities created by recent inf all/mergers (e.g.. lAdami et al 



2004), although density 



20051) are likely to influ- 



ence the pressure maps of the Coma cluster that is dynamically young. Obviously, there are 
other potential driving sources including AGN and subcluster mergers, ram pressure strip- 
ping, etc. It will be interesting to see how turbulence driven by each process adds together 
when nonlinear effects as well as magnetic fields are considered. 
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Fig. 1. — (a) Profiles of the radial orbit (tUorb.r), circular orbit {uJorh,c), and Brunt- Vaisala 
(cubv) frequencies in our cluster model, (b) X-y plane snapshot at t = 6 Gyr of the per- 
turbed density relative to the initial value in logarithmic scale for a circular-orbit galaxy, (c) 
Averaged input rates of the kinetic energy due to a single galaxy as functions of its mass. 
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Fig. 2. — Temporal evolution of the kinetic energy inside the cooling radius for adiabatic 
models that differ in the number and mass of the galaxies. Solid curves show the results 
from runs with 256^ zones. Dotted lines (for models D and E) from 128^-zone runs are 
within 8% of the corrresponding solid curves. Dashed and dot-dashed lines draw Ek = 
8.2 X 10^^M^^^(t/ Gyr) erg for Mn = 1 and 5, respectively. 
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Fig. 3. — (a) Three-dimensional velocity dispersions cr^ of the ICM turbulence driven 
by galaxy motions as functions of NgMf^. Open circles, with errorbars representing the 
standard deviations in the temporal fluctuations of a^, denote the density- weighted val- 
ues, while the non-weighted ones arc given by filled circles. The solid line corresponds to 
ay = 2.2 km s'^iVg^^Mn. (b) Spherically binned power spectra of (solid line) compressive 
and (dotted line) shearing parts of the ICM velocity in models C and E at t = 6 Gyr. 



